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We examine one important (and overlooked in all previous investigations) aspect of well - known 
crossing diabatic potentials or Landau - Zener (LZ) problem. We derive the semiclassical quantiza- 
tion rules for the crossing diabatic potentials with localized initial and localized or delocalized final 
states, in the intermediate energy region, when all four adiabatic states are coupled and should be 
taken into account. We found all needed connection matrices and present the following new analyti- 
cal results: (i) in the tunneling region the splittings of vibrational levels are represented as a product 
of the splitting in the lower adiabatic potential and the non-trivial function depending on the Massey 
parameter; (ii) in the over-barrier region we find specific resonances between the levels in the lower 
and in the upper adiabatic potentials and in that condition independent quantizations rules are not 
correct; (iii) for the delocalized final states (decay lower adiabatic potential) we describe quasi - 
stationary states and calculate the decay rate as a function of the adiabatic coupling; (iv) for the 
intermediate energy regions we calculate the energy level quantization, which can be brought into 
a compact form by using either adiabatic or diabatic basis set (in contrast to the previous results 
found in the Landau diabatic basis). Applications of the results may concern the various systems, 
e.g., molecules undergoing conversion of electronic states, radiationless transitions, or isomerization 
reactions. 

PACS numbers: 31.50.Gh, 05.45.-a, 72.10.-d 



I. INTRODUCTION. 

There has been great progress in the theory of crossing potentials during the last seven decades (see e.g. the 
references in both research and textbook literature [1] - [9]). Surprisingly, seemingly a simple but a basic question 
how well known semiclassical quantization rules should be modified for this particular situation (crossing diabatic 
potentials with bound (i.e. localized) initial and localized or delocalized final states) to the best of our knowledge 
are still unanswered, (at least a complete and unifying description of the quantization for a general case is still not 
available and a number of other questions remain to be clarified). 

Partially it is related to the fact that unfortunately experimental data in this field are still scarce and not very 
accurate. However, the situation is now changing. Experimental techniques (e.g., the increasing precision of experi- 
mental tests in the femtosecond laser pulse range enables to excite well defined molecular states and to study their 
evolution in time using the second probing laser beam [5], [10], [11]) have progressed to the point where molecular 
tunneling splitting dependence on energy can be measured in well controlled conditions with a very high accuracy. It 
would therefore seem appropriate at this time to take a fresh look at the theory, which is the purpose of the present 
article. Note also that these questions are not only of interest in their own intellectual right. Recent experimental 
and theoretical advances [12], [13] in particular are beginning to yield a coherent understanding of several phenomena 
that, far from requiring minor corrections to the standard adiabatic treatment of the problem. Physically such kind of 
situations can occur as a result of non-adiabatic interactions of different electronic states forming in crossing one-well 
diabatic vibrational potentials. Adiabatic coupling removes diabatic level crossing, and the diabatic levels are replaced 
by the adiabatic ones (see Fig.l illustrating this phenomenon). In the case of a large adiabatic splitting (see precise 
criteria below) one can restrict oneself to the only lower adiabatic potential (symmetric or asymmetric double-well, 
or decay potential for the systems under consideration) and neglect any influence of the upper adiabatic potential 
(parabolic one- well for our case). However, in a general case of arbitrary adiabatic splittings, intra- well and inter- wells 
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dynamics depends on the both adiabatic potentials (i.e. on tunneling and adiabatic splitting). Having in mind appli- 
cations, the studies of these questions may concern the various molecular systems undergoing so-called conversion of 
electronic states, isomerization reactions, or radiationless transitions arising from "intersystem" crossings of potential 
energy surfaces in molecular spectroscopy and chemical dynamics, or inelastic atomic collisions. 

It is worth noting that there exists a huge literature devoted to different approaches have been made by other 
authors to the problem of crossing diabatic potentials (see e.g. [2], [3], [6] - [10], [14] - [20]) but some important 
differences to our work should be noted. First, suffice it to say, that the problem how diabatic potentials crossing 
modifies the adiabatic potentials (occurring as a result of this crossing) quantization rules, has not been investigated 
at all. One of the reason is that for many say standard rigid molecules with quite large adiabatic splitting of energy 
levels, one may safely neglect any influence of the upper adiabatic potential (i.e. to use the standard quantization 
rules). However, nowadays the increasing precision of experimental tests of molecular tunneling splitting and decay 
(and besides investigations and synthesis of more and more new non-rigid molecules), makes the study of this problem 
relevant and actual. Second, is a methodic note. All previous approaches were based on the general semiclassical 
WKB formalism. The crucial point to treat quantization for crossing diabatic potentials is how to compute the 
contribution coming from the contour around a complex turning point. The accuracy of the WKB method can be 
improved considerably, [14] - [17] (more recent references on so-called Laplace contour integration can be found also 
in [2], [13]) by the appropriate choice of the integration path around the turning point. This method ascending to 
Landau [1] appears to be quite accurate for the tunneling and over-barrier regions, (however, even in this case there 
are some non- negligible corrections found in the papers [6], [18], [19]), but in the intermediate energy region (where 
there are relevant contributions from all four quantum states occurring at the crossing diabatic potentials) the method 
becomes completely non-adequatc. Besides the choice of these additional special trajectories (which one has to include 
to improve the accuracy of the WKB method, and along which the semiclassical motion is described by the Weber 
functions) depends on the detail form of the potential far from the top, and therefore for each particular case the 
non-universal procedure should be perform from the very beginning. 

The essential simplification of the procedure can be achieved using the standard WKB semiclassical approximation 
in the momentum space representation was also proposed in the literature [2], [6], [15] - [19]. The method works well 
to compute the LZ transition probability, however, the application of this approach to the level quantization problem 
is difficult to realize. Indeed the problem requires to know the eigenfunctions in the coordinate space, and one can 
not bluntly use the Fourier transform of the functions found in the momentum space, since the WKB method gives us 
only asymptotics of the eigenfunctions. These say drawbacks of the WKB-like methods did not allow to study level 
quantization for crossing diabatic potentials in the previous investigations, and, besides, we believe we are the first 
to explicitely addresses the question on the behavior in the intermediate energy region. In all previous publications 
this region was considered as a very narrow and insignificant one, or in the best case the results were obtained by a 
simple interpolation from the tunneling (with monotonic decay of the transition probability) to the over-barrier (with 
oscillating behavior) regions. 

Recently we have shown [21] that semiclassical solutions of many eigenvalue problems can be considerably simplified 
by including into the consideration second order turning points. The fact is that one second order turning point replaces 
two close linear turning points. Moreover it turns out (see below) that connection matrices which link on the complex 
plane the solutions to the Schrodinger equation in the vicinity of the the crossing point with the asymptotic solutions 
far from this point, can be calculated from the solutions of the Weber equation. Increasing and decreasing solutions 
in the classically forbidden region around a second order turning point, are characterized by the action which has 
a minimum along a certain trajectory, we will refer in what follows as the instanton trajectory. The same kind of 
the Weber equation can be formulated to calculate the connection matrices in the vicinity of a saddle point (or a 
maximum) of the potential, but besides we should also relate increasing and decreasing solutions at the crossing point 
(see also our recent publications on LZ crossing phenomena [22]). 

Our aim is to construct semiclassical wave functions. To do it we use connection matrix methodology which can be 
applied to any semiclassical approximation, but details of the method depends on the order of the turning points. For 
the second order turning points which are minima of a potential, the whole procedure is equivalent to the traditional 
instanton approach, and the imaginary time (i.e. after Wick rotation) instanton trajectories correspond to the periodic 
orbits between the turning points, and the connection matrices in this case (see below and appendix A to the paper) 
are real-valued ones. It is not the case for the second order turning point which is the potential maximum. The 
complex- valued connection matrix links two regions of infinite motion. Formally one might refer the corresponding 
wave functions also as instanton ones applying twice the Wick rotations. LZ crossing points are combinations of 
two second order turning points with two different Stokes constants [23] corresponding to one minimum and one 
maximum of the potential. In the tunneling region there exist periodic orbits for two solutions while two others 
correspond to unrestricted (infinite) motions. As above we will call these wave functions as instanton ones since they 
are Weber functions (like in the traditional instanton method) but with complex-valued arguments occurring as a 
result of complex coordinate frame rotation. 
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Thus our approach in this paper is based on the minimization of the functional of classical action in the upside-down 
potential, so-called instanton type approach, which represents the most important area of the configurational space 
where the semiclassical wave functions are localized ( see [4], [13] - [22]). The whole analysis can be brought into a 
more elegant form by introducing connection matrices which link on the complex plane the semiclassical solutions to 
the Schrodinger equation for the model potential of the problem under study and the exact solutions of the so-called 
comparison equation which is valid near the crossing point, where one can approximate the potential by linear or 
parabolic ones. The explicit calculations of the connection matrices are rather involved since one should treat the 
four fundamental solutions to the left and to the right regions with respect to turning or crossing points. Therefore 
the connection matrices, we are looking for, are 4x4 matrices. Although the generalization for our case of the 
known already 2x2 connection matrices (see e.g., [23]) is straightforward, it deserves some precaution as it implies 
quite different procedures for the energy (more accurately for E/j were 7 ^> 1 is the semiclassical parameter, see 
below) smaller (the tunneling region), larger (the over-barrier region), or of the order (the intermediate region) of 
the potential barrier energy. Within the framework of the connection matrix approach we present a full and unified 
description of II? (which is very often can be quite reasonable approximation for real systems) level quantization 
problem for diabatic potentials crossing. 

The remainder of our paper has the following structure. Section II contains basic methodical details and equations 
necessary for our investigation. Except for a mathematical trick that eliminates a large amount of tedious algebra 
and allows us to construct a regular method for calculating higher order perturbative corrections, the section contains 
already known results. New physical results are collected in sections IV, V and VI, and partially in section III, where 
we calculate all needed connection matrices, which provide a very efficient method of finding semiclassical solutions 
to the Schrodinger equation in potentials having several turning points. The knowledge of the connection matrices is 
important and significant not only in itself but also for developing a good analytical approximation, and standardized 
numerical procedures. In Section IV we find the quantization rules for the tunneling and over-barrier energy regions. 
Section V is devoted to the intermediate energy region where all four states occurring at the diabatic potentials crossing 
should be regarded on the same footing. Different particular cases, depending on the ratio of the model parameters 
are also examined in this section. In section VI we investigate the linear coupling of the LZ system to harmonic 
phonons and find that it renormalizes of the parameters entering the initial diabatic potentials crossing problems 
considered in the previous sections IV - V. The last section VII contains summary and discussion. In two appendices 
to the paper we collect some more specialized technical material required for the calculations connection matrices in 
different energy windows (appendix A), and to reduce fourth order Schrodinger equation to two independent second 
order Weber equations (appendix B). Those readers who are not very interested in mathematical derivations can skip 
these appendices finding all essential physical results in the main text of the paper. 



II. FOURTH ORDER COMPARISON EQUATION FOR THE CROSSING POINT. 



To move further on smoothly let us describe first our strategy. First we should define all notations and relevant 
points of the diabatic potentials crossing problem. We depicted the typical situation in the vicinity of the diabatic 
potentials crossing point in Fig. 1. The diabatic potentials (1,2) are shown by thin solid lines, the adiabatic potentials 
(3,4) by bold solid lines. Besides we have introduced in the picture the adiabatic coupling energy U12, the crossing 
point energy U# , and E is the characteristic zero-point oscillations energy in the parabolic barrier approximated the 
lower adiabatic potential near its top. 

As a model for diabatic potentials in this paper we choose two parabola Ul, U r with a symmetrical crossing in the 
point x = 0. To be specific: let us consider two types of the diabatic potential crossing depicted in the Fig. 2. The 
corresponding adiabatic potentials are, respectively, the double well or decay lower potential, and the one-well upper 
adiabatic potential. At arbitrary values of the parameter U12 to find eigenstates and eigenfunctions for our model 
potential we should solve the coupled Schrodinger equations 



ld 2 0. 



2 dx 2 



+ ~/ 2 (U L (x) - E)Q L = 7 2 u 12 e R ; - 



i d 2 e R 

2 dx 2 



+ 1 2 {U R {x)-E)Q R = 1 2 U 12 Q L . 



(2.1) 



Here 7 ^ 1 is the semiclassical parameter which is determined by the ratio of the characteristic potential scale over 
the zero oscillation energy (i.e. 7 = mQ aQ/fi, where m is a mass of a particle, a is a characteristic length of 
the problem, e.g. the tunneling distance, f^o is a characteristic frequency, e.g. the oscillation frequency around the 
potential minimum). 

These equations (2.1) can be written as one fourth order equation, e.g. for Ql 



d 4 e L 

dx* 



-2 1 2 (U L (x) + U R (x)-2E)^± 



^ 2 du L de L 

dx dx 



-4 7 4 



(U L -E)(U R -E)-U- 



12 



1 d 2 U L 
2 7 2 dx 2 



Ql =0(2.2) 
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In what follows we use ilo and ao to set corresponding dimensionless scales, e.g dimensionless energy e — E/jClo, 
u l/r = Uj^/ji/^hCto), U12 = 2Ui2 / (~fh£l Q ) (we introduce factor 2 in u 12 for ease of writing following below equations), 
coordinate X = x/ao, and we put K = 1 (except where explicitely stated to the contrary and the occurrences of h are 
necessary for understanding). 

Luckily the equation (2.2) admits semiclassical solutions by Fedoryuk method [27] - [29] since the coefficients at the 
n-th order derivatives proportional to 7~™, and therefore so small that higher order derivatives of the prefactor (in 
the semiclassical form the wave function can be always presented as the prefactor times the exponent) can be safely 
neglected in finding of asymptotic solutions. Fedoryuk method makes possible to find asymptotic solutions to the 
ordinary differential equations of the following form 



,(«) 



ra-l 
fc=0 



7*(*)0 (fc) = o, 



(2.3) 



where we designated j/ fe ) = d k y/dX k , and the coefficients at the derivatives fk(X) are arbitrary functions of X. Note 
that Eq. (2.2) for has this Fedoryuk form. By the standard semiclassical substitution y = A exp(-yW (X)) (2.3) 
can be reduced to the set of equations combining the terms proportional to 7™, 7™" 1 • • • , which for 7 >• 1 can be 
represented in the form of generalized so-called Hamilton - Jacobi and transport equations, respectively 



n—k 



F(\)=\ n + J2fk(X)\ k =0, 



fc=0 



and 



dF dA 
IXdX 



1 d 2 F dX 
2~dJ?Jx 



A = 0. 



(2.4) 



(2.5) 



where A = -~/dW/dX. 

Noting that in the vicinity of the crossing point X = the parabolic diabatic potentials can be replaced by the 
linear ones counted from the barrier top C/ # 



u L/R (X)=u*±fX, 



(2.6) 



(as above w # = C/ # / (7MI0)) , and eventually the equation (2.2) can be presented into a more compact and simple 
form 



dX 4 



L r, 2 

— 27 a 



d*e L 2 d@ L 
dx 2 1 J dX 



- 277^ + 7> 2 - fx 2 uf 2 ]e L = , 



(2.7) 



where in our dimensionless units a = 2(u# — e). 
The roots of the characteristic polynomial for (2.7) 

F(X, X) = \ i - 2 7 2 K + u R - 2e)X 2 - 4 7 2 ^A + 4 7 4 



(u L - e)(u R - e) - u\ 2 - 



1 d 2 u L 



2 7 2 dX 2 



= 0, 



(2. 



or in the equivalent form 

F(A) = A 4 — 2a 7 2 A 2 - 2 7 2 /A + ^(a 2 u\ 2 f 2 X 2 ) , (2.9) 
determine independent solutions to (2.7). Solving the equation F{\) = perturbatively over 7 _1 <C 1 we find 

Xj = A^° + A) , (2.10) 

where 

ni/2 



A = ±7 



a± 



l 12 



J2 X , 



(2.11) 



and 



± 



2^1 



f 2 X 2 



(2.12) 
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we find finally the four asymptotic solutions of (2.7) 



{ Vj } = {QX, e;,e+e:} = ( — ) exp 



dF 



-1/2 



r \ 3 (x')dx> 

Jo 



(2.13) 



The subscripts in (2.13) corresponds to the upper or lower adiabatic levels, and the superscripts are referred to the 
sign of the action. 

As it was mentioned above in the vicinity of the crossing point one can replace (2.2) by (2.7), and by the substitution 



we can find the equation for k 



where 



e L = exp(K li2 X)$^ 



k - a 7 V + - 7 X 2 = -k 4 5 2 (1 + 26) + R{k, 6) , 



(2.14) 



(2.15) 



R(k, S) = (2k 6 )- 1 (1 - 3<5)(1 + 5)~ 3 (l -Q- ^1 - 2Q 2 )) ; Q = 85 2 (1 + 5) , 



7 2 f 1 

and 

Ki, 2 = ±lV^ (l - y) • (2-17) 

It can be proved that the 4-th order equation (2.7) in variables (2.14) is reduced to the equation with constant 
coefficients in front of all derivatives and with the free term in the form of a quadratic over X function. In the case 
when the exponent in (2.14) is a solution to the equation (2.15)thc transformed equation is reduced to the Weber 
equations upon neglecting anharmonic terms like X 2 dF/dX, X 3 F(X), and X 4 F(X). We presented all details of this 
reduction in appendix B to the paper. Thus the equation (2.7) is reduced to two independent Weber equations with 
the known fundamental solutions [34] 



{6 L } = jexp(±7v^0I>-„ (j^pj ' xj > exp(± 7> /5A")I>-i-„ (j^pj 1 xj j - ( 2 - 18 ) 



where v = "fu\ 2 / '(4/ 'y/a) is so-called Massey parameter. The corrections to the indices of the parabolic cylinder 
functions D and to the arguments of these functions can be found from (2.15). 

Presented above the leading terms of these solutions corresponding to the tunneling case, i.e. (here we use a 
dimensional energy E) 

E<(U* - U 12 ) (2.19) 

(in our dimensionless units it is a > U12), where the characteristic fourth order polynomial (2.8) can be reduced to 
the second order one (i.e. two pairs of roots are nearly degenerated), are known in the literature (see e.g. [1] - [4]) 
but the Fedoryuk method we used, gives us also in the tunneling region the higher order over the parameter 5 (2.16) 
corrections. 

In the tunneling region (2.19) one can expand the roots of (2.15) in terms of the parameter 5 (2.16). Using 
the substitution (2.14) to transform (2.2) we can find easily that at the conditions (2.15), (2.16) the coefficients at 
the fourth and at the third order derivatives in the transformed fourth order differential equation for $ are small 
(proportional to S and to VS respectively) and thus this fourth order equation can be rewritten as two second order 
Weber equations with the solutions 

D pll ,v(0X), 

where 



4 4 " V a J V 4 
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The same manner in the over-barrier energy region i.e. 

E>(U# + U 12 ) (2.21) 

(again as in (2.19) we have used dimensional units, and in dimensionless form (2.21) reads as —a > U12), when the 
energy is larger than the upper adiabatic potential minimum, the roots of the equation (2.8) arc complex - conjugated 
and having the same structure as presented above (see also (B6) in the appendix B) for the tunneling region with the 
roots k given 




Besides in (B6) 



'V2 V 

and with S playing the role of the small parameter in this region 

~ 6 =^W (2 - 23) 

Again as above for the tunneling region, the coefficients at the higher order derivatives are small, and therefore, the 
function $ (2.14) satisfies the Weber equation with the fundamental solutions 



£>p<i,2) (/3(i,2)A) 



where 



5 S 
p (1) = -1 -i- + iO, p {2) = i- + iO, (2.24) 



/ 2f2\l/4/ ?2 \ / 2 £2 

■ ' • l 2 ^) I T ,A=exp(-i37r/4)fy 




(D is defined as the Massey parameter entering (2.18) with a — ► \a\, i.e. v = {iu\ 2 / (Af yj\a\)). Like it was for the 
tunneling region (2.20), the leading terms of the expansion (2.24) coincide with the well - known results, but from 
(2.24) we are able to compute the corrections to the main terms. 

The analogous task for the intermediate energy region, i.e. (in dimensional units) 

{U* + U 12 ) >E>(U#- C/12) , (2.25) 

is much more tricky. Our results will be presented in section V, but a few comments are necessary here. In the 
problem we have three dimensionless parameters characterizing the energy (a), the level coupling (1*12), and the 
potential (/) , and besides for the ease of semiclassical estimations we keep also the semiclassical parameter 7 > 1. 
Note also that these parameters are not independent ones, and the relation u\ 2 = 2/ 2 (which we will be useful in our 
further consideration) should be satisfied. In terms of these parameters within the intermediate energy region (2.25), 
we have the subregion, S', \a\ < Zj^ 1 , and U12 < 27 _1 , and the intermediate subregion, 5", where 7^/^12/2 3> 1. In 
the section V we calculate the connection matrices for the both subregions, and details of the reduction procedures, 
which are different in S' (where the comparison equations are reduced to two decoupled Airy equations) and S" 
(where these comparisons equations are Weber ones) are described in the appendix B. 



III. CONNECTION MATRICES. 



The purpose of this section is to briefly indicate the main steps in the derivation connection matrices. The matching 
points we must find to quantize the energy levels depend essentially on the energy window under consideration (2.19), 
(2.21), (2.25). The tunneling region is placed in the lower adiabatic potential. In the WKB method in this case 
starting from the crossing point (X = 0) one has to investigate the classically forbidden region in the lower adiabatic 
potential barrier (see Fig. 3a and the corresponding figure caption for all notations) . The solutions can be found easily 
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in the vicinity to the crossing point but to derive the quantization rules one should know also the solutions quite far 
from the crossing point. To do it explicitely in the WKB method we should match the two exponentially decreasing 
and two exponentially increasing solutions in the barrier with the oscillation solutions in the wells. Technically the 
matching should be performed asymptotically, i.e. at small \X\ but for large enough ^/t"|X|. To do it one has to 
calculate all needed connection matrices (namely at the crossing point, and at the linear and second order turning 
points, and the shift matrices from the crossing point to the turning points in the classically forbidden region and 
between the turning points in the classically accessible region). Within the instanton type method the trajectory goes 
through only the classically forbidden region (see Fig. 3b), and matching should be performed only at two second 
order turning points. In the over-barrier energy region (see Fig. 4), the matching is performed by using the crossing 
point connection matrix U' c , the shift matrix L connecting the crossing point and the linear turning points at the 
upper adiabatic potential, and the shift matrices L L / R belonging to the diabatic potentials. In this case all matching 
solutions are oscillating ones. Finally for the intermediate energy region no real- valued turning points for the upper 
states (see Fig. 5) and the matching between two oscillating and two exponentially varying solutions is determined 
by the connection matrix (see Fig. 5) linking the linear imaginary turning points of the adiabatic potentials. 

To recast the analysis into a compact form it is convenient to formulate the general procedure for calculating of the 
connection matrices for an arbitrary combinations of the first and of the second order turning points. After that the 
procedure can be applied to any particular problem under investigation. To do it, one has to extend the known for 
linear turning points procedure [23] . For a generic semiclassical equation 

— +1 2 q {z)^(z) = Q, (3.1) 

in the limit 7 >• 1 the Stokes and anti-Stokcs lines are determined by the following conditions, respectively 

RW(z) = 0, (3.2) 

and 



$SW(z) = 0, (3.3) 

where the action 

W(z) = J ^{z)dz. (3.4) 

The lower integration limit in (3.4) is not relevant because we are interested in semiclassical solutions for large 
\z\. These Stokes and anti-Stokes lines separate the complex plane z into the sectors. On the anti-Stokes lines the 
increasing and decreasing solutions become equal, and the type of the solutions is interchanged upon crossing of the 
anti-Stokes lines. The Stokes lines are bisectors between neighboring anti-Stokes lines. After the crossing with the 
Stokes line, one should add to the coefficient at the decreasing solution the coefficient at the increasing solution times 
so-called Stokes constant. The latter one is occurred as a result of going around the turning point and depends on 
the turning point type. 

To find the connection matrices for the tunneling region we have to establish the correspondence between the 
solutions of the fourth order differential equation (2.2) and those for the localized in the left (L) and in the right (R) 
wells states. In the case a 3> f\X\ for the diabatic potentials, the action can be computed starting from the both 
wells (R and L) 

^Wl ~ 7 ^oL + k a X + ^X 2 , jW* R 7 VF * K - k X + ^X 2 , (3.5) 

where ko = jy/a is imaginary momentum ( U# sets the energy corresponding to the diabatic potentials crossing 
point), and W ' are the actions computed from an arbitrary distant point in the L or in the R wells, respectively 
to the point X = 0. From the other hand in the adiabatic potentials — u# ± \Ju\ 2 + f 2 X 2 the corresponding 
actions can be represented 

B 2 

1 W ± ~ jW± = k X ± '—X 2 sign{X) . (3.6) 

Explicitely comparing the semiclassical wave functions in the both representations (adiabatic and diabatic ones) it is 
easy to see that the adiabatic functions in the potential U~ coincide with the diabatic functions for localized L and 
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R states at X < and X > respectively. The adiabatic functions for the upper potential U + correspond to the 
tails of the diabatic wave functions localized in the opposite wells. Therefore in the level crossing region the L/R 
diabatic functions are transformed into the R/L functions, and the interaction entangles the diabatic states with the 
same sign of koX. Thus we have only four non-zero amplitudes of the following transitions 



Recalling that 



1 W ± = 7 / (a ± 



1/2 



k Q x±P-x*± v -(i 



In v) 



(3.7) 



(3i 



we come to the conclusion that the quantum solutions (2.20), valid in the vicinity of the level crossing point asymp- 
totically, match smoothly increasing and decreasing solutions, and it leads to the Landau description [1] of the level 
crossing transitions. To illustrate the presented above analytical results we show schematically in Fig. 6 the matching 
of the asymptotic (Fedoryuk) solutions (2.13) for the crossing diabatic potentials with localized initial and final states 
via the Weber functions in the tunneling region. We use the symmetric basis constructed from the functions (2.13) 
(see detail description in the Fig. 6 caption). 

In the tunneling region (2.19) for every well (L or R) there exist increasing and decaying exponentially real-valued 
solutions to the Schrodinger equation. The solutions are matched at the crossing point, therefore they are linked 
by the real-valued 4x4 connection matrix which should have two 2x2 blocks linking the increasing (decreasing) 
diabatic solution in the L-well with the decreasing (increasing) diabatic solution in the i?-well, in the agreement with 
the standard Landau scheme of the tunneling transitions [1]. Omitting a large amount of tedious algebra we can 
represent the connection matrix linking the "asymptotic" (i.e. in the left/right (L, R) wells and for the upper/lower 
(+, — ) adiabatic potentials) solutions in the tunneling energy region in the following form 



M C (+) L^ C) M C ( - 




U c 



F c M { c +) L { 1 ] M [ c 




(3.9) 



Here U c is the 4x4 connection matrix at the crossing point, which in the tunneling region has the following form 

(3.10) 



p 

(sin 2 (7w))/p — cos(7r^) 
cos(nv) p 
cos(7w) 



— cos(7w) 



(sin 2 (7ri/))/p 



where we designated 



P : 



2t t exp(-2x) 
I» 

r(+) 



(3.11) 



and x — — (1/2) (y — (1/2)) In v. The matrices M^' and AfJ •* are the 2x2 connection matrices at the 

corresponding turning points, which are determined by the phase shifts at these points 



Mi") = 



1 -i \ 
-(i/2) (1/2) 



(3.12) 



and Mc + ^ is the matrix Hcrmitian conjugated to (3.12). The L^J R and F c matrices are called shift matrices, and 
those are related to the variations of the coefficients of increasing and decaying semiclassical solutions in the regions 
between the turning points (F c is the shift matrix when one moves from the crossing to the turning point in classically 

forbidden region, and L^} R are the shift matrices in the classically accessible regions). Explicitely we get 



( exp(- 7 W£/2) 

exp( 7 W|/2) 



(3.13) 



Here Wg is the action in the lower adiabatic potential barrier. Finally the structure of the shift matrices L^j R is 



He) _( exp(i 7 VF£ /fl ) 

L LIR - I 





c M-nWi /R ) 



(3.14) 
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where W£, R is the action calculated by the integration between the turning points. We present explicit expressions 
for the total connection matrix in the appendix A. 

The same manner can be treated the over-barrier region (2.21). In this case the crossing point is in the classically 
accessible region for the both potentials. The fundamental diabatic solutions can be represented as the waves prop- 
agating in the opposite directions, and the complex-valued connection matrix has as it was for the tunneling region 
2x2 block structure, where the blocks link the waves in the L and in the R wells propagating in the same direction. 
Specifically the corresponding connection matrix at the crossing point U' c 



U' 



sexp(— i(f>) — exp(— ttv) 

sexp(i^) — cxp(— irv) 

exp(— ttv) sexp(—i<p) 

exp(— nv) sexp(i(p) 



(3.15) 



where we denoted s = \J\ — cxp(— 2nu), = argT(— iv) + 3(2%), and 

X = -(i/2)((7r/4) + - In i/)) + (1/4) (tti/ + Ini/) , 



(3.16) 



should be multiplied by two blocks: the block from the left gives the contribution at the turning point and includes 
the shift matrix to the crossing point in L and in R wells of the lower adiabatic potential; the right block is related to 
the turning point and to the shift matrix to the crossing point in the upper one-well adiabatic potential. Thus finally 
in the over-barrier region we get 



M c (+) i (c) 



* - u - \u' 
M(+)L J c 







(3.17) 



Here we used the same notations as it was above for the tunneling region, and besides the matrices are transposed 

with respect to the matrices Mj^ given in (3.12), and the new shift matrix L is 



exp(-i 7 T^*/2) 

exp(ijW*/2) 



(3.18) 



(remind that W* is the action in the upper adiabatic potential). Combining altogether (3.12) - (3.18) one can trivially 
find the full connection matrix for the over-barrier energy region (2.21). We present the explicit form of the matrix 
in the appendix A. 

More tricky task is to calculate the connection matrix in the intermediate energy region (2.25). Following the same 
line as above we first present the general structure of the connection matrix in the intermediate energy region 



o M ( r } 



u" 



m| +) 



(*t\ 



(3.19) 



These matrices M±^ have been introduced in our paper [21] for the imaginary turning points characterizing the both 
adiabatic potentials in the region |a| < U12, and they read 



1 

(i/2)cxp(- 7 Wf) 



where Wf~ are so-called Euclidian actions in the reversed upper and lower adiabatic potentials, which can be estimated 
as 

1 ~ 7 ' 

where 



91,2 = 



1U\2\/U12 ± a 



10 



and all other matrices entering (3.19) are defined already. Since M±^ turn into the unit matrices at a > U12 and 
a < U\2, the connection matrix (3.19) matches continuously into the corresponding matrices (3.9), (3.17) for the 
tunneling and over-barrier regions. 

The connection matrix in the intermediate energy subregion 5"' can be calculated using Weber function asymptotic 
expansion for large complex indices (B15), which are the solutions to the comparison equations (2.15) in the inter- 
mediate energy subregion S" . These 4 roots are arranged clock- wise and counter clock- wise on the radius 7-^/^12/2 
circle around the crossing point. The following combinations of the comparison equations match the semiclassical 
solutions (2.13) 



+ ej <-» e j + 63 <-» e+ er; + et + e 4 



e: 



Combining together the asymptotic expansions for these combinations, we find at the crossing point, the matrix U" 



is 



U'J 



(V2n/T(q*))cM-2x(q*)) _ 

(T(q)/V2^) exp(2 X («))(l - cxp(-2^ 2 ) cos 2 (^i)) 
exp (— 27r<j , 2 )cos(7rg 1 ) 

exp ( — 2nq2 ) cos (irqi ) 



(3.20) 



exp ( — 2nq2 ) cos (irqi ) 

— exp(— 2nq2) cos(7rgi) 
( v / 2^/r(g))exp(2 X (<z)) _ 

(IV)/v^F)exp(2x(g*))(l -cxp(-27rg 2 )cos 2 ( 7 rg 1 )) 



where as above 



4/ 

and, besides, we introduce the following abridged notations 



q = qi + iq2 ; qi,2 = 77 ; <? =qi~iq 2 , (3.21) 



and analogously 



\ = \ 1 - '\2 : 2\ 1 = - ( (71 - - ] In Y f ^ • (3.22) 



2X2 = qi - qi In \q\ - <p ( qi - - ) , (3.23) 



where ip is defined by (BIO). Now the full connection matrix in the both intermediate energy subregions can be found 
easily simply collecting the given above expressions, and the explicit form for the connection matrix is presented in 
the appendix A. Note that the intermediate energy region connection matrix (3.20) has the same block structure as 
the connection matrices in the tunneling and in the over-barrier regions. This is a consequence of the fact that in a 
neighborhood the diabatic potential crossing point, only the Weber functions with equal indices can be hybridized. 
At \a\ = U12 the connection matrix (3.20) turns into the connection matrices (3.10) for the tunneling region and into 
(3.15) for the over-barrier energy region, and it enables us to construct semiclassical solutions for any arbitrary energy 
window. Note, however, that in the intermediate energy region the Massey parameter is replaced by the complex 
index q. In the section V we will present another derivation of the connection matrix (3.20), and will discuss specific 
relations between the adiabatic and diabatic states in the intermediate energy region. 



IV. QUANTIZATION RULES. 

In the tunneling energy region, one has only real-valued eigenfunctions, since in the both wells there are only 
the localized states. In this energy window the connection matrix linking the "asymptotic" (i.e. in the left/right 
(L, R) wells and for the upper/lower (+, — ) adiabatic potentials) solutions is represented in the form (3.9) and 
is given by Exps. (Al) - (A7). Within the WKB method we should match the two exponentially decreasing and 
two exponentially increasing solutions in the barrier with the oscillation solutions in the wells, thus it requires the 
knowledge of the connection matrices at the crossing point and at the linear turning points, and the shift matrices 
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from the crossing point to the turning points in the classically forbidden region and as well between the turning points 
in the classically accessible region (Fig. 3a). Within the instanton type method the trajectory goes through only the 
classically forbidden region (see Fig. 3b), and to perform the matching one should know also the connection matrices 
for the second order turning points. 

The same manner can be treated the over-barrier region with the corresponding connection matrix (3.17), and (A9). 
Evident modifications of the given above expressions for the tunneling and over-barrier regions should be performed 
to treat the intermediate energy windows. Indeed in this case one has to take into account also the contributions 
from the imaginary turning points. The procedure is reduced to replacement of oscillating factors by exponentially 
decaying ones (see details in the next section). Finally for the intermediate energy region no real- valued turning points 
for the upper states (see Fig. 5) and the matching between two oscillating and two exponentially varying solutions is 
determined by the connection matrix (3.19) linking the linear turning points of the adiabatic potentials (see Fig. 5). 



A. Two diabatic parabolic potentials crossing. 

Now (collecting the explicit expressions for all needed connection matrices from the appendix A to the paper) we 
are in the position to derive the quantization rules, which can be formulated as a condition that the amplitudes of 
exponentially increasing at X > 0, and X < 0, respectively, solutions must be vanished. Taking into account 

that Wl = Wft (the actions in the corresponding wells of the lower adiabatic potential) and using the connection 
matrix relating the fundamental solutions of the Weber equation, we can formulate the corresponding quantization 
rule for the tunneling region in terms of the matrix elements defined by (3.9) 

^22^33 - m 23 m32 = , (4.1) 

where are corresponding matrix elements from (3.9). 

Putting all together we can find from (3.9) - (3.10), and (4.1) the quantization rule for this case 

tan( 7 W*) = ±|exp( 7 ^), (4.2) 

where Wg is the action in the barrier formed in the lower adiabatic potential, and p = U\\ is the corresponding matrix 
element of the connection matrix (3.10). 

Only the factor l/p varying from to 1 in the diabatic and in the adiabatic limits, respectively, makes this 
quantization condition (4.2) different from the well known [1] quantization rule for the symmetric double- well potential. 
Correspondingly, the tunneling splitting at finite values of the so-called Massey parameter v can be represented as a 
product 

A n = A° n p(v), (4.3) 
of the tunneling splitting A° in the adiabatic potential and the factor 

/2tt 



P(" ) = f(^~ (1/2) exp(-i/) , (4.4) 

associated with the transition amplitudes between the diabatic potentials in the crossing region. 

It is particularly instructive to consider (4.2) as the standard [1] Bohr-Sommerfcld quantization rule, where in the 
r.h.s. the both, geometrical ip n and the tunneling \n phases are included additively. In the adiabatic limit when 
p(v) — > 1, we find that ip n ^> and (4.2) is reduced to the quantization of the symmetric double-well potential. In 
the diabatic limit ip n = — Xn and the geometric phase compensates the tunneling one. The physical argument leading 
to this compensation may be easily rationalized as follows. Indeed, at the reflection in the crossing point X = 0, the 
trajectories in the classically forbidden energy region are the same as those for the tunneling region but with a phase 

Shift 7T. 

We focus now on the quantization rules for the over-barrier energy region. Closely following the consideration 
performed above for the tunneling region, and replacing the connection matrix (3.10) by the corresponding matrix 
for the over-barrier region (see section II) (and making some other self-evident replacements) we end up after some 
algebra with the quantization rule 

(1 - exp(-27w)) cos(2 7 W^ - <p) cos( 7 IT* + <p) + exp(-2nv) cos 2 ( jW[ + ^— — ] = , (4.5) 
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where W* is the action in the well formed by the upper adiabatic potential, and <j> = argT(— iv) + 3(2%) is determined 
according to (3.16). From the Eq. (4.5) follows that the eigenstates are determined by the parameter 



B = 



exp(— 2-kv) 
1 — exp(— 2ttv) 



(4.6) 



In the diabatic limit v — > 0, and therefore, B — > l/(27r^) in (4.5) the main contribution is due to the second term, 
and it leads to a splitting of degenerate levels in the diabatic potentials. Moreover since 



7 



W* L 



VP 
~2~ 



7r n + 



- ± 2f sin 
2 



W* 



(4.7) 



the splitting increases when the Massey parameter v increases, and it is an oscillating function of the interaction U\2- 
In the adiabatic limit, when v — > oo, — > 0, and, therefore, from (4.6) B ~ cxp(— 27w), the main contribution 
to (4.5) comes from the first term which determines the quantization rule for the upper one-well potential and for 
the lower double-well potential in the over-barrier energy region, and in this limit the parameter B plays a role of 
the tunneling transition matrix element. For B smaller than nearest level spacings for the lower and for the upper 
potentials, one can find from (4.5) two sets of quantization rules leading to two sets of independent energy levels 



jW* = 7T ( Til 



1 



; 2jW[ 



77 n 2 



(4.8) 



Since the eigenstate energy level displacements depend on the adiabatic coupling U12 the resonances can occur at 
certain values of this parameter, where the independent quantization rules (4.8) are not correct any more. The widths 
of these resonances are proportional to exp(— 2m>) and therefore are strongly diminished upon the Massey parameter 
v increase. This behavior is easily understood, since in the limit the wave functions of the excited states for the 
lower potential are delocalized, and their amplitudes in the localization regions for the low-energy states of the upper 
potential, are very small. 



B. Bound initial and decay final states: the diabatic potentials (l + X) 2 /2 and (1/2) —X crossing. 

The second instructive example treats the one-well and linear diabatic potentials crossing. It leads to the lower 
adiabatic decay potential and to the upper one- well adiabatic potential. The quantization rules in this case correspond 
to the vanishing amplitudes for the exponentially increasing solutions when X — > — 00, and besides one has to require 
that no waves propagating from the region of infinite motion, i.e. at X > 1/2. Performing the same as above procedure 
we find that in the tunneling energy region, the eigenstates are the roots of the following equation 

tan( 7 T4^) = -i-i- exp(2 7 W* ) , (4.9) 

with the same as above notation. To proceed further it is convenient to introduce the complex action to describe the 
quasi - stationary states 

« = *(f- ! i)' ( " 0) 

where evidently fi = OWl/OE does depend on E. From (4.10) the real and imaginary parts of the quantized 
eigenstates are 

E n = n (n + ; T n = p 2 (^ cxp(-2 7 T4^) . (4.11) 

This relation (4.11) describes the non-adiabatic tunneling decay of the quasi-stationary states of the lower adiabatic 
potential. The same as we already got for the two parabolic potentials crossing (4.3), here the tunneling and the 
adiabatic factors are entering decay rate multiplicatively. Since the decay rate is proportional to the square of the 
tunneling matrix element, T n oc v 2 (y) as it should be. 
In the over-barrier energy region the quantization rule is 

(1 - cx P (-2tw) cxpH(7VF£ - </>)] cos( 7 VF* + <p) + (4.12) 

exp(-27w) exp(-ijW*/2) cos ( 7 VF£ + ) = , 
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and the actions depend on the energy E as 



"fW'l = ir— ; 7VF* = 7r —7 



n (u* + u 12 ) e_ 



(4.13) 



where f2 and £~2i are independent frequencies of the diabatic and the upper adiabatic potentials. 
In the diabatic limit the decay rate is proportional to the Massey parameter v and has a form 



r„~-^cos 2 ( 7 w* + 4>), 



(4.14) 



and in the opposite, adiabatic, limit the decay rate is 



r„ ~ n exp(-27ri/)(l - sin(2 7 IT£ - <t>)) . 



(4.15) 



In the both limits the decay rate is the oscillating function of Uu. 

We illustrate the dependence Y(U\2) for the crossing diabatic potentials U\ = (1 + X) 2 /2 and U2 = (1/2) — X in 
the Fig. 7. Note that while the tunneling decay rate of the low-energy states is increased monotonically with the 
Massey parameter v, the decay rate of the highly excited states goes to zero in the both (diabatic and adiabatic) 
limits. Besides there are certain characteristic values of U\2 when the r.h.s. of (4.14) or (4.15) equal to zero and 
therefore r„ = 0. This seemingly paradoxical and contradicting to conventional wisdom result can be rationalized as 
follows. For the case under consideration (one well upper adiabatic and decay lower adiabatic potentials) there are 
always energy levels blocked by the upper adiabatic potential. This resonance phenomenom manifests itself wave- 
like particle properties omnipresent in quantum mechanics. For the system under consideration the upper adiabatic 
potential is equivalent to a resonator with a set of well - defined modes (resonances) with high quality factors. An 
important feature (in distinction to a conventional resonators where these modes occupy more or less homogeneously 
the whole phase space) is that the resonance modes are localized in its own effective cavity whose position is given by 
the conditions r„ = (4.14) or (4.15). 

Quite similar one can study the more general example, describing two non-symmetric diabatic potentials crossing 
at X = point: 



In a certain sense it is the generic case, and when the parameter b entering the potential (4.16) is varied from 1 to 
00, we recover the two particular examples considered above, and come from two identical parabolic potentials to the 
case one- well and linear diabatic potentials crossing. This kind of the potential was investigated recently by two of the 
authors (V.B. and E.K) [30] aiming to study crossover behavior from coherent to incoherent tunneling upon increase 
of the parameter b, the larger is this parameter b, the larger will be the density of final states. The criterion for 
coherent-incoherent crossover behavior found in [30] based on comparison of the transition matrix elements and the 
inter level spacings in the final state. The analogous criterion should hold for the level crossing problem, however in 
the latter case the tunneling transition matrix elements has to be multiplied by the small adiabatic factor. Therefore 
the coherent - incoherent tunneling crossover region moves to the more dense density of final states, and the larger U12 
is the smaller will be the region for incoherent tunneling. Quite different situation occurs for highly excited states. 
In the diabatic limit, the transition matrix element is increased with the Massey parameter v, and therefore at a 
given b value, the system moves to more incoherent behavior. In the adiabatic limit, the transition matrix element is 
exponentially small, and coherence of the inter-well transitions should be restored. However, since the matrix elements 
are oscillating functions of U12 for the intermediate range of this coupling (U12) coherent - incoherent tunneling rates 
are also non-monotonically varying functions. 



More difficult task is to derive the quantization rule in the intermediate energy region, where all four roots of 
the characteristic equation contribute into the solutions. One has to use the connection matrix (3.19) computed for 
this region (see details in section III and appendix A). It has two 2x2 blocks structure, the same as the connection 
matrices for the tunneling and over-barrier regions. Here, we present another derivation of the same connection matrix 
using the adiabatic representation. It offers a deeper insight into the mathematical structure of the problem, and 
besides provides physically relevant relations between the adiabatic and diabatic states in the intermediate energy 
region. The very possibility to use the both representation is stipulated by the fact (we have mentioned already 



-(l + Xf-u 2 = -{X*-2bX + b). 



(4.16) 



V. 



INTERMEDIATE ENERGY REGION. 
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in section III) that the semiclassical eigen functions in the intermediate energy region can be represented as linear 
combinations cither diabatic or adiabatic functions (this adiabatic - diabatic transformation has been discussed for 
quantum coherence phenomena in [31], see also [2]). 

Since the adiabatic potentials have two second order turning points (the minimum of the upper, and the maximum 
of the lower adiabatic potentials) the blocks of the connection matrix in the intermediate energy region (where 
now unlike the matrices (3.9), (3.17) describing the transitions between the diabatic states, the connection matrix 
corresponds to the transitions between the adiabatic states, and non-adiabatic perturbations induce the transitions 
[21]), are characterized by the parameters qi i2 analogous to gi ; 2 from (3.21) entering (3.20). Respectively, for the real 
- valued blocks 

7\/2wi2 , . , . 

9i = — jj — ( M i2 + a )' I 5 - 1 ) 

and the complex - valued blocks ( associated with the maximum of the lower adiabatic potential) 



g2 = 4/ (ui2 - a) . (5.2) 

We can now reap the fruits of the previous subsection consideration efforts. First, let us note that from the relations 
(2.20) and (2.24) one can see that when the energy approaches to the top of the barrier, the exponents and pW 
of the parabolic cylinder functions are increased and thus, more and more deviated from the value prescribed by the 
Massey parameter v. Second, increasing of f}^) upon |a| decreasing, decreases the values of \X\ where the asymptotic 
smooth matching of the solutions should be performed. For 5 — > these \X\ values are located deeply in the classically 
forbidden region, where the potentials are close to the diabatic potentials, while for 5 > 1/4, these coordinates \X\ 
are of the order of the quantum zero-point oscillation amplitudes, and therefore to find the solution in this region, we 
have to use the adiabatic representation. 

Although as it is shown in the appendix B, the intermediate region for the both subrcgions, 5"', at S < 1/4, and S', 
at 5 > 1/4 can be investigated on equal footing in the frame work of the comparison equations (i.e. at the diabatic 
basis) it is instructive to study the problem also in the adiabatic representation, what is the purpose of this section. 
As a sub-product of this consideration we get also the justification of the comparison equation approach. In the 
adiabatic basis the intermediate subregions S' and S" should be studied separately. Two simple observations give us 
a conjecture how to treat the problem in the intermediate energy region. First of all the energetical " window" for the 
intermediate subregion 5"', where 6 < 1/4, and \a\ < U12, in terms of the dimensional energy scale is determined by 
the rectangle around the crossing point 

U12 < 2C/r 2 ; \U* -E\< , (5.3) 

where we define U{ 2 = (l/2)(h 2 F 2 /to) 1 / 3 . By the other words the characteristic interaction energy at the intermediate 
region boundaries does not depend on ?7i2. Analogously the intermediate subregion S' is restricted by the lines 

U12 > 2C/r 2 ; \U* -E\< U 12 . (5.4) 

The positions of the linear turning points \X*\ corresponding to the energies U# ±Ul 2 do depend on the ratio U\ 2 /Ul 2 . 
These points are located inside or outside of the interval [— 7 -1 ^ 2 ; +7~ 1//2 ] at Ui 2 /Ui 2 < 1 and at U\ 2 /Ul 2 > 1, 
respectively. Accordingly for the both cases the matching conditions in the intermediate energy region are different. 
In the former case for the asymptotic matching region the potentials can be reasonably approximated by parabola, 
and therefore we should work with the Weber equations, and for the latter case the matching are performed in the 
region where the potentials are linear ones, thus the equations are reduced to the Airy ones. 

Let us discuss first the intermediate energy subregion S" , where q\ and q 2 are large, and therefore the Massey 
parameter, i.e. the indices of the Weber functions are also large. The arguments of the Weber functions are oc X^fj, 
and their asymptotic expansions determine the interval where the matching should be done (5.4). In what follows 
we will closely follow the method we borrowed from Olver paper [32] (for the asymptotic expansions of the Weber 
functions with large indices, see also his monograph [33]), which is in fact an expansion over small parameters l/|g,| 
(where \(ji\ are the exponents (5.1), (5.2)) of the fundamental Weber solutions, and it leads to the following asymptotic 
solution at X > 

#+(x) ~ r+ 1/2 (a + r+)-« i ex P (- 7 xy + ) , *:(x) ~ y: 1/2 {x + r_) l « 2 ex P (i 7 Ay_) , (5.5) 



where Y± = \Ju\ 2 ± a 2 + f 2 X 2 . Using the known relation between the fundamental solutions of the Weber equation 
[34] 



£> M (z) = rxpi ;-//:■/.),,•::.' • —^7 exp '~ l ~~r~ ) ' ' 
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we can find two other (complimentary to (5.5) solutions 

= Y~ 1/2 (- sin(7rgi)(X + cxp(- 7 Xr + ) + e xp(-2 X i) f ^y|^y (X + Y+)« cxp( 7 Xr+)^ , (5.6) 

and 

*±(X) = (5.7) 



Y+ 1/2 (-*exp(-7rg 2 )(X + Y_)^ cxp^XY^) + exp(-2 X2 ) r((1/ ^ ..^ (X + Y_)^ expH 7 XY_)j . 

In the case of weak level coupling, i.e., for the intermediate energy subregion S" , the adiabatic potentials everywhere 
(except a small neighborhood of the level crossing point) can be linearized, i.e. represented as a ± f\X\, and the 
asymptotic solutions are reduced to a linear combination of the following functions 

$± cx (f\X\)-V 2 cxp(±e+ + s^nX) , <f>± « (/l^l)- 1 ^ cxp(±C- - signX) , £± = ^j(.f\X\ ± a) 3/2 , (5.8) 

and these functions are smoothly matched with semiclassical solutions (see details in the appendix B). As a result we 
can calculate finally the connection matrix U'J in the intermediate energy region in the adiabatic basis 



(^/r(-zq 2 ))exp(-2 X ( J q 2 )) 

(r(«7i)/V2^)cxp(2 X (<?i))sin 2 (^ 1 ) 

cos(7r<7i) 

iexp(— Trq 2 ) 



-2CXp(-7T<7 2 ) 

— cos(— qi) 
(v^F/r(g 1 ))exp(-2 X (gi)) 

2(T(—iq2)/V2n) cxp(-2ix(g 2 )) exp(-7rq 2 ) cosh 2 (7rg 2 ) 

where the function \ i s defined in (3.22), (3.23). We see that the connection matrix in the adiabatic basis, unlike 
(3.20) defined in the diabatic basis, does not provide continuous transformation into the connection matrices for the 
tunneling and over-barrier energy regions ((3.10) and (3.15) correspondingly). This apparent inconsistency is due to 
disregarding of adiabatic level interactions which become relevant in the intermediate energy region. However there 
is a simple remedy to ensure the continuous over all energy windows matching of the connection matrices. One has to 
rotate the complex plane q over the angle (BIO). Thus luckily (as it is often the case in semiclassical approaches) 
we can safely reduce the problem quite accurate to the Weber or Airy equations in the both intermediate energy 
subregions, using respectively the perturbation theory with respect to the diabatic or adiabatic states. The found 
adiabatic connection matrix could be used on the same footing as the diabatic connection matrix (3.20), e.g., to derive 
the quantization rule, which for the intermediate energy window can be written in the simple and compact form as 

cos(2 7 T^£) = - cxp(-7T(j 2 ) . (5.9) 

It is useful to illustrate the essence of the given above general result by simples (but yet non trivial) examples. 
First, let us consider two identical parabolic potentials with their minima at X = ±1 and with the coupling which 
does not depend on X. Since the symmetry, the solutions of the Hamiltonian can be represented as symmetric and 
antisymmetric combinations of the localized functions 

* ± = 4s(*i ±$ «)- ( 5 - 10 ) 
V2 

The functions are orthogonal, and, besides, two sets of the functions , \&jT ), and (^q > ^e) (where the subscripts 
and e stand for the ground and for the first excited states respectively) correspond to the two possible kinds of level 
crossings. 

In Fig. 8 we depicted schematically the dependence of the level positions on the coupling [7 12 . In the energy region 
E < U* + U\2 where only there exist the discrete levels of the lower adiabatic potentials, there are the pairs of the 
alternating parity levels (^+ , *q ), and (^q , ^~). The tunneling splittings are increased monotonically since the 
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Massey parameter v is increased, and the barrier is decreased with [7 12 . The same level and parity classification is 
remained correct for the energy region above the barrier of the lower adiabatic potential where the spectrum becomes 
almost equidistant one. However, in the over-barrier region, the resonances are occurred between the levels of the 
same parity, and this sequence of the odd and of the even levels is broken, and level displacements are not monotonic 
functions of Uw Some of the levels of different parities can be mutually crossed. For the upper adiabatic potential the 
level sequence is opposite to this for the lower adiabatic potential. The intermediate subregion S" limits are shown 
by two dashed lines. The boundaries between the intermediate subregion S' and the tunneling and the over-barrier 
regions are shown by the dotted-dashed lines outgoing from the corners of the subregion S" rectangle, and these lines 
coincide with energetic displacements of the top and of the bottom of the adiabatic potentials. Note also that we 
checked the results of our semiclassical approach and found remarkably good agreement with the numerical quantum 
diagonalization. Shown in the Fig. 8 level displacements versus U\i coincide (with the error not exceeding 10% for 
the full range of variation of f/12, including the both intermediate energy subregions) with the results of the numerical 
diagonalization in the basis of harmonic oscillator functions of the initial Hamiltonian (2.1 for two diabatic crossing 
potentials (1 ± X) 2 /2. 



We have considered semiclassical quantization of bound and quasi-stationary states beyond the adiabatic approxi- 
mation but for ID case only. Of course the energetic profile of any real system is characterized by a multidimensional 
surface. However, it is often possible to identify a reaction coordinate, such that the energy barrier between initial 
and final states is minimized along this specific direction, and, therefore, effectively one can treat the system under 
consideration as II?, regarding all other degrees of freedom as a bath of harmonic oscillators. In this section we 
investigate the simplest multidimensional Hamiltonian describing the non-adiabatic transitions, namely the 2x2 ma- 
trix potential for the X variable (or what is the same two ID diabatic potentials crossing considered in the previous 
sections) and the set of "transverse" harmonic oscillators {Yfc} coupled with the reaction coordinate X 



Here V\ (X) is the bare (in a general case anharmonic) ID potential, cok is the eigenfrequency of the transverse oscillator 
k, the function F(X) describes how the only strongly fluctuating coordinate X is coupled to thermal bath of transverse 
oscillators, and Ck are corresponding coupling constants. This kind (6.1) of the multidimensional potential has been 
studied in the literature (see e.g. [4]), and some efforts were made to find a feasible approximation to treat the 
potential within the semiclassical approach. In this section we legitimate the method proposed in [4] focusing on 
the LZ problem in the tunneling region. Similar consideration can be easily generalized for the over-barrier and 
intermediate regions. 

The equations of classical motion (in imaginary time) for the transverse coordinates have the following form 



VI. COUPLING TO A THERMAL RESERVOIR. 



,2 




(6.1) 



fc 



fe 




k 



(6.2) 



where III is the integral transformation 




(6.3) 



It can be expand in the following high- and low-frequency limits 




(6.4) 



where 




(6.5) 



At the high frequency limit (6.4) is reduced to the trajectory equation but with the renormalized potential corre- 
sponding to the following X-dependent effective mass 
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where 




(6.7) 




(6.8) 



In the p4 approximation for the spectral density of oscillators the neglected in (6.7) terms proportional to p 6 ( an d 
the last term in (6.8)) are small. The physical message of the calculation performed in this section is that the 
renormalization of the effective mass leads to slowing down of the motion and it is equivalent to say that the Massey 
parameter is renormalized 



where X c is the crossing point. Of course the coupling will also change the action along the extremal action trajectory 
(this effect has been discussed in the literature, see e.g. [4]). The specific for the LZ problem new phenomcnom is 
renormalization of the Massey parameter (6.9) which controls main features of the behavior for any system undergoing 
level crossing. 



In conclusion we end up stressing again the main point of our methodology. We have shown that the comparison 
equations for the 4-th order differential Landau-Zener equations in the coordinate space can be represented as two 
decoupled Weber equations. The indices and the arguments of the corresponding Weber functions defined by the roots 
of the characteristic equation (2.15) for the complex wave vector k, and \k\ 3> 1 in the semiclassical approximation. 
In the frame work of our method the diabatic potential crossing points are treated as two second order turning points 
characterizing by different Stokes constants [23]. The accuracy of the method depends on anharmonic terms, which 
are not taken into account in the comparison equations, but which are small in the semiclassical approach over small 
parameters S, S, or 5i nt respectively, in the tunneling, over-barrier and intermediate subregion S" energy windows. In 
the subregion S', 5i nt is not a small parameter. However, since the asymptotically smooth matching is performed at 
small \X\ < 7~ 1,/2 , anharmonic corrections to the comparison equations can be safely neglected for this subregion as 



We have presented detailed semiclassical analysis of crossing diabatic potentials problem. We examine one important 
(and overlooked in all previous investigations) aspect of well - known energy level quantization problem for crossing 
diabatic potentials. We derive the semiclassical quantization rules for the particular situation of crossing diabatic 
potentials with localized initial and localized or dclocalizcd final states, in the intermediate energy region, when all 
four adiabatic states are coupled and should be taken into account. In fact it exhausts all cases practically relevant 
for spectroscopy of non-rigid molecules (i.e. with more than one stable configuration). 

We use the connection matrix methodology which presents a simple and standardized description of any semi- 
classical approximation, which offers therefore a deeper insight into the mathematical and physical structure of the 
approximation. We found that in the tunneling region the tunneling splitting is represented as a product of the split- 
ting in the adiabatic potential and the non-trivial function p{v) (we calculated analytically) depending on the Massey 
parameter, i.e. on the energy and the slopes of the diabatic potentials in the crossing region. In the over-barrier 
region we found specific resonances between the levels in the lower and in the upper adiabatic potentials and in that 
condition one may not use independent quantization rules. New results have emanated from our consideration of the 
intermediate energy region. For this energy region we calculated the energy level quantization, using adiabatic basis. 

We have presented in this paper all details of the LZ problem for two electronic states using the connection matrix 
approach for the LZ problem in the coordinate space, the approach which turned out very efficient for this class of 
problems, and are important in many areas of pure and applied sciences. Even though only model potentials are 
investigated here, our approach is quite general and has potential applicability for various systems in physics and 
chemistry, and the results can be tested by their experimental consequences for many examples of molecular systems 
undergoing conversion of electronic states, non-radiative transitions, or isomerization reactions, and not only. The 
results of the LZ-problem investigations are very relevant for slow atomic or molecular collisions , [35], [2], where the 
interaction of diabatic potentials induces transitions between initial and final electronic states. However, since the 



v — > v* = v\J 'm*(X c ) 



(6.9) 



VII. CONCLUSION. 



well. 
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interaction is essential only near the crossing point, one can compute the transition probability, linearizing the both 
diabatic potentials (see our consideration in section II). The same approximation works quite well for the so-called 
pre-dissociation phenomena. 

However (in contrast to the atomic and molecular collision problems) there are fundamental problems of chemical 
physics and molecular spectroscopy where one may not restrict oneself to the only transition probability calculations, 
but should know the complete cigcnvalucs/cigcnfunctions solution. It is the case for example if we are interested in 
the calculation of vibrational - tunneling spectra of non-rigid molecules, or reactive complexes with more than one 
stable configuration. The lowest multi - well potential of such systems is formed from one well diabatic potentials 
crossing corresponding to each stable configuration. Apart from the lowest potential, the upper adiabatic potential 
with its minimum above the maximum of the lowest potential should be also taken into account for these situations 
(see Fig. 1). In the most of the calculations of tunneling splittings in the ground and low excited vibrational states 
the coupling to the upper potential are neglected, what is certainly correct only for strong enough adiabatic coupling. 
Evidently it is not the case for the levels close to the adiabatic barrier top, and especially in the upper potential well. 
The quantization of these levels play noticeable in the spectroscopy of non-rigid molecules, and the same situation 
takes place for systems undergoing the Jahn - Teller effect, where the interference of the diabatic states occurs in this 
energy region [3] . 

One more example for the application of our results is molecular radiationless transitions within excited electronic 
states. Typically for this situation the decay potential is formed owing to crossing of bound and unbound diabatic 
potentials. Since the radiationless transitions are followed by luminescence and chemical reaction phenomena (see e.g. 
[36] - [39]) one should know the complex eigenvalues of the quasistationary states prepared by optical pumping. 

Let us also stress that in real systems the characteristic values of the coupling between the diabatic states can vary 
within the very wide range from several eV for the electronic states of the same symmetry to zero (for the states 
with different spins) . To treat all these cases one should know the solution of the diabatic potentials crossing problem 
described in our paper for the corresponding wide range of the Massey parameter from v — to v 3> 1. 



The research described in this publication was made possible in part by RFFR Grants. One of us (E.K.) is 
indebted to INTAS Grant (under No. 01-0105) for partial support. The authors thank also their manuscript referees 
for insightful criticism. 



Putting all given in the section II expressions (3.9) - (3.14) together we can recapitulate the matrix elements rriij 
of the full connection matrix in the tunneling region 
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APPENDIX A 



mn = -. exp(— 7W5) cos(7W£) cos^W^) exp^Wg) sin^W^) sin(jW^) ; 




(Al) 



(A2) 



(A3) 



(A4) 



(A5) 
(A6) 



(A7) 



sin 2 (itv) 

m 3 3 = P ; W44 = ; = m 43 = 

V 



(A8) 
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For the over-barrier region the full connection matrix could be given in a more compact form. Using (3.17), (3.15), 
(3.18), (3.12) from the main body of the paper we get the following matrix 



(s/2) cos(jW L r - 4>) s siii^Wlr ~ cj>) - cxp(-7w) sm(^W Rit ) - exp(-7w)/2 cos(^W R *) 

— ssin(-fWLR — 4>) 2s cos^Wlr — </>) — 2 cxp(— nu) cos^Wr*) exp(— ttv) sm(qWR^ 

— cxp(— -kv) sin(7VFi *) 2 exp(— irv) cos^V^l ») 2s cos^W* + <f>) — s sin^W* + (f>) 

exp(-7w)/2cos(7Wi,*) exp(— 7rz/) sin^Wz, *) s sin(-fW* + <p) (s/2) cos^W* + </)) 



(A9) 



where W LR = W* L + W R , and W L/R> , = W* L/R + W*/2. 



APPENDIX B 



The efficency of the standard instanton approach [25], [26] (see also [4], [22]) is based on a successful choice of 
the comparison equation near second order turning points, where asymptotically smooth matching of semiclassical 
solutions to the solutions of this equation should be performed. It is known for example [21] that for anharmonic 
potentials the Weber equation provides such a very successful choice since in the matching region anharmonic cor- 
rections are still small. The aim of this appendix is to show that the analogous situation holds for crossing diabatic 
potentials points, where two Weber equations can be successfully used as the comparison equations to the fourth order 
Landau-Zencr equation (2.2). The arguments and the indices of the fundamental solutions to these Weber comparison 
equations are determined by the roots of the corresponding characteristic equations (see below and the main body of 
the paper). 

To prove the statement let us first substitute (2.14) into the equation (2.7). We get 

L> 4 $ + 4kL> 3 $ + (6k 2 - 2a 7 2 )D 2 $ + 4(k 3 - a 7 2 K - i 7 2 /)L>$+ (Bl) 



[k a - 2a~/ 2 K 2 - 2 1 2 f K + 7 4 (a 2 - w 2 2 - / 2 X 2 )]$ = 



where D n = d n /dX n . The equation (Bl) can be formally derived by simple manipulations (two sequential differenti- 
ations and summations) from the following second order equation 



where the coefficients are 



£> 2 $ + (a + aiX + a 2 X 2 )<S> = , 



7 2 f 

a = K 2 - a 7 2 -J-L(l + S); ai = J 2 f5 ; a 2 = --/ 2 fn8 , 



(B2) 



(B3) 



where n should be found from the characteristic equation (2.15), and 8 is given by (2.16). 
The fundamental solutions to (B2) read as 



± 



7 4 r 1/4 



X - 



2k 



where 



P - 



+ 



7 4 / 2 



-1/2 



a 



4a 2 



(B4) 



(B5) 



In the tunneling (2.19) and over-barrier (2.21) regions of energies these 4 solutions (2 solutions of (B4) for two largest 
modulus roots of the characteristic equation (2.15)) can be separated into two independent pairs. In the tunneling 
region the two largest modulus roots of (2.15) are (two other roots are small and do not satisfy semiclassical approach) 



k = Kq I 1 ± 



8 2 



2 2k 2 q - aY 



; 



_L 

V2 



+ \J~a 2 



'12 



1/2 



(B6) 



Putting (B6) into (B5) we find (neglecting S 2 terms, i.e., for k = k ) 4 fundamental solutions to the comparison 
equation in the form (2.18). Thus from the given above expressions and (2.16), (2.17) from the main text we conclude 
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that the solutions &l/r (2-18) can be expanded over our small parameter 6, and due to the condition (2.16) anharmonic 
corrections to the Weber functions (B4) are small (by other words the parameter 5 determines the accuracy of our 
approximation) . Indeed the anharmonic terms neglected in the Weber comparison equations are of the order of 5 (it 
is an upper estimation at X = a/f, i.e. at the boundaries of the intermediate energy region), thus the corrections are 
small according to (2.16). The same kind of analysis can be performed in the over-barrier region (2.21), where one 
finds two imaginary largest modulus roots of the characteristic equation. The roots are given by (B6) with kq and 
the small parameter 5 defined according to (2.22), (2.23). 

One simple observation helps to perform the same analysis for the intermediate energy region (2.25). Indeed, since 
the differences between the solutions to the characteristic equations for A (2.9) and for k (2.15) determine the accuracy 
of our approach, let us compare the solutions. The roots of (2.9) at X = 

Ai,2 — ±7\/a + U12 ; A 3 , 4 — ±7\/a - uli (B7) 

are moved upon the variation of a in the intermediate energy region from the real to imaginary coordinate axis. 
Analogously the roots of (2.15) 

1/2 / N 1/2 



«i,2 - ±-^= [a + \Ja 2 - uj 2 J ; k 3 . a ~ ±-^= [a - \j a 2 - uj 2 J (B8) 

are moved along the real and imaginary axis in the tunneling and in the over-barrier regions respectively. 

We conclude from the (B7) and (B8) that in the tunneling and in the over-barrier energy regions there is one-per-one 
correspondence between the roots A of (2.9) and k of (2.15). Just this correspondence allows us to match smoothly the 
semiclassical solutions to the Schrodinger equation and the Weber functions found as the solutions to the comparison 
equations. It is not the case in the intermediate energy region where two roots of (2.9) are real and two are imaginary 
ones having the same modulus, i.e. moving upon a variation along a circle with the radius 7-^/^12/2. In this case 
the semiclassical solutions can be presented as certain linear combinations of the comparison equation solutions. We 
have found these combinations in the adiabatic basis in the section V. In this appendix we show how to solve the 
same problem in the diabatic basis, and it reveals more clearly and explicitely an estimate of the omitted terms in the 
equation and the areas where the solutions become wrong and where the matching procedure is carried out. Indeed 
the roots of (2.15) in the intermediate energy region (2.25) are 

«i,2 - ±7^/^- exp(i<£) ; n 3A ~ iry^^ exp(-itp) , (B9) 

where 



tan<^ = W — -. (BIO) 

V "12 + a v ' 

Correspondingly to these roots (BIO) the arguments and the indices of the Weber functions (B4), (B5) read as 

z 1 = z 2 = 2K int y / ~5~~ t exp(-i<p/2)(X + (2K int )~ 1 exp(-iip) ; (Bll) 



zz = z±=- 2n mt y/S^t exp(i(p/2)(X + (2K int ) 1 cxp(i(p) , 



and 



Pi = Pi ~ 1 



4S: 



int 



■ exp(-i<p)(l + Sf nt exp(-2i<p)) ; 



(B12) 



Pa = P3 ~ 1 = -1 - tt — exp(i<p)(l + S mt exp(2i^)) , 

where K int = 7(ui 2 /2) 1/2 , and S int = (7 2 /)/(4«L)- 

Using known due to Olver ([32], [33]) asymptotics of the Weber functions we are in the position to compare 
the semiclassical functions with the solutions to the comparison equations. The former functions determine by the 
exponential factor 

7/ 2 



F± (X) = 7 V^±^X + X 3 , (B13) 

12mi2a/mi2 ± a 
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while the exponential factors entering corresponding asymptotics of the Weber functions are 

7/ 2 



iq, 2 (A) = 7 V^J±^(1 + S int )X ± Kf nt 8f nt exp(-2ip)X 



12zti 2 V u i2 ± a 



1 ± <>int 

«12 



X 6 



(B14) 



Let us consider now the intermediate subrcgion 

5", M < (f/l) 2/3 , and u 12 < 2/7, (see (5.3)), where (2.16) does not 
hold. Luckily, however, the asymptotically smooth matching is performed at small \X\ < 7 -1 / 2 , where the comparison 
equation (B2), and, therefore, the characteristic equation (2.15) arc valid (although, S is not a small parameter). In 
this subregion we have to take into consideration the term R(k, S) in (2.15). At a = 0, and iti 2 = 0, the characteristic 
equation has one double degenerate root k = 0, or correspondingly in (B2), a 2 = 0. Thus the comparison equations 
are reduced to two decoupled Airy equations. Using known Olver asymptotics for the Weber functions with large 
arguments and indices [32], [33] 



D p (z) oc exp 



1/2 



dz 



we can find asymptotics to the solutions of (B2) 

$0 oc cxp i J \/ clq + a\X + aiX 2 dx 



(B15) 



(B16) 



valid at arbitrary values of the parameters (a 2 = including). This relation (B16) provides asymptotically smooth 
matching of the semiclassical solutions with the Weber functions in the intermediate subregion 5"' (where k is of the 
order of 7 3> 1, and with the Airy solutions in the subregion 5", when k ~ ^7. 

This consideration provides the justification of our approach described in the main body of the paper. As it is seen 
from (B13), and (B14) at small a the accuracy of the asymptotically smooth matching of the semiclassical solutions 
with the Weber functions is of the order of 5i ntl and close to the energetic boundaries (2.25) of the intermediate 
region, anharmonic corrections (A 3 ) are increased. Thus we conclude that the matching for this case (2.25) can be 
performed either in the adiabatic basis (as it has been done in the section V) or in the diabatic basis as we have 
shown in this appendix. The simplest way to prove the equivalence of the both representation is to transform into 
exponential forms the factors like (X + Y + ) qi etc, entering the solutions (5.6), (5.7), found in the section V. In the 
both methods the accuracy is of the order of Si nt , and the connection matrices presented in the appendix A do not 
depend on the basis. 
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Figure Captions. 

Fig. 1 

Potentials in the vicinity of the diabatic potentials crossing point 

The diabatic potentials (thin lines, 1,2), the adiabatic potentials (bold lines, 3,4) by bold solid lines, the adiabatic 
coupling energy U\2, and E is the characteristic zero-point oscillations energy in the parabolic barrier approximated 
the lower adiabatic potential near its top. The tunneling energy E region is shown by a broken line. 

Fig. 2 

The diabatic level crossing phenomena 

(a) bound initial and final states; 

(b) bound initial and decay final states. 
Fig. 3 

Connection matrices for the tunneling energy region: 

(a) in the WKB approach to the lower states, where are the connection matrices for the linear turning points, 
and U c for the crossing point; the shift matrices are depicted as arrows, in the classically accessible regions Ll and 
Lr, and in the classically forbidden region Fl and Fr (for the upper states no real- valued turning points); 

(b) in the instanton type method one has two connection matrices M L j R for the second order turning points and 
shift matrices Fl and Fr in the classically forbidden region. 

Fig. 4 

Connection matrices for the over-barrier energy region. The shift matrices from the crossing point to the inner 
turning points are designated by L (all other notations are the same as in the Fig. 3). 
Fig. 5 

Connection matrices for the intermediate energy region (like in the tunneling region no real valued turning points 
for the upper states). 
Fig. 6 

The matching of the asymptotic solutions in the tunneling region for the diabatic levels crossing shown in Fig. la: 
1 - the function <S>+(X)V2^/T(l + v); 
2- the function <3>£(X); 

3 - the function <Z>r~(X) y/2n/T(l + v); 

4 - the function ^(X); 

V - the function exp(fc X)D_i_ I/ (/3X); 
2' - the function ex.p(k X)D-- L _ l/ (-(3X); 
3' - the function exp(— koX)D-i- u (f3X); 
4' - the function exp(-k o X)D_ 1 _ v (-0X). 
Fig.7 

T„ versus U\2 for the quasi stationary states at the diabatic potentials (1 + X) 2 /2 and (1/2) — X crossing; (a) 1 
- 4 are the level energies 0.042 , 0.125 , 0.208 , and 0.292 for the lower adiabatic potential; (b) V — 3' are the level 
energies 0.625 ; 0.708 ; 0792 for the upper adiabatic potential. 

Fig.8 

Level displacements versus U\i for two diabatic crossing potentials (1 ±X) 2 /2. Dashed lines show the intermediate 
energy region (the subregion S" is between the dashed lines, while the subregions 5" are confined to the left pockets 
between the dashed and dotted-dashed lines); dotted - dashed lines show also displacements for the top and for the 
bottom of the adiabatic potentials, k, n, and n' are quantum numbers for the diabatic, and lower and upper adiabatic 
potentials. Note that shown in the figure level displacements coincide with the error not exceeding 10% with the 
results of the numerical diagonalization in the basis of harmonic oscillator functions. 
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